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Abstract 

We consider nonlinear parabolic equations involving fractional diffusion of the form 
dtu + (—A)®$(m) = 0, with 0 < s < 1, and solve an open problem concerning the 
existence of solutions for very singular nonlinearities <f> in power form, precisely = 

gy-(n+i) fQj. gome 0 < n < 1. We also include the logarithmic diffusion equation 
dtU + (—A)'* log(M) = 0, which appears as the case n = 0. We consider the Cauchy 
problem with nonnegative and integrable data uo(x) in one space dimension, since the 
same problem in higher dimensions admits no nontrivial solutions according to recent 
results of the author and collaborators. The limit solutions we construct are unique, 
conserve mass, and are in fact maximal solutions of the problem. We also construct 
self-similar solutions of Barenblatt type, that are used as a cornerstone in the existence 
theory, and we prove that they are asymptotic attractors (as t —>■ oo) of the solutions 
with general integrable data. A new comparison principle is introduced. 
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1 Introduction 

In this paper we consider a class of nonlinear parabolic equations involving fractional dif¬ 
fusion of the form 

(1.1) (9iUA)^$(n) = 0 . 

The symbol (-A)^ denotes the fractional Laplacian operator with 0 < s < 1, i. e., the 
nonlocal operator defined by 

(1.2) {-Ayv{x) = c{N,s)p.v. [ \N+ 2 s ^y^ Vx€R^, 

Jw.N \x - 

acting on the whole Euclidean space at least for functions in the Schwartz class S. The 
formula is valid for all dimensions N > 1. The constant c{N,s) is given in the literature 
but it is not needed in what follows and p.v. means principal value of the integral. 

The existence and properties of solutions for this type of equations with fractional diffusion 
has been studied by the author and collaborators for nonlinearities that are positive and 
increasing for u > 0, in particular when <h(u) = with m > 0, cf. [TSl HU 1201 El 1431142]. 
This includes singular cases for 0 < m < 1 since then <h'(?x) = mu^~^ —>■ oo as n —)• 0. 

Here, we are interested in very singular nonlinearities, more precisely, when <1> : M_|_ ^ M 
is a monotone increasing function of u with a singularity in u = 0 such that <I>(0-|-) = — oo. 
Consequently, nonnegative data and solutions are considered. The standard cases we have 
in mind are 4>„(u) = — 1/u"' for some n > 0, or <I'o(^^) = log(u). They correspond to 
4>'(u) ~ with n -|- 1 > 1, thus the denomination very singular introduced in the 

literature for this type of equations with standard Laplacian, cf. [39] • We will keep this 
tradition for equations with a fractional Laplacian. These very singular diffusion equations 
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are also described in the literature as very fast diffusion, superfast diffusion, or ultra-fast 
diffusion, cf. e. g. [2511311139]. 

For such equations the existence of solutions is not at all obvious. Thus, we have proved 
in a recent paper with Bonforte and Segatti [7] that when the space dimension is iV > 2 
and we try to solve the Cauchy problem in the whole space with integrable initial data, 
then there exist no nontrivial solutions, even if we accept local-in-time solutions defined for 
a short time interval, 0 < t < T0 The same happens for the problem posed in a bounded 
domain with zero Dirichlet data. 

The purpose of this paper is to prove that there is a range of existence of solutions for very 
singular equations of the form with ‘h = if the space dimension is 1. We will also 
prove that the solutions have the good properties of the non-singular range of parameters 
$(u) = with m > 0. A very crude explanation of the existence result is as follows: by 
becoming strictly positive for t > 0, the solutions avoid the singularity in a way that suffices 
to grant first nontrivial existence, and then the rest of the properties. 

Before stating the results, let us point out that the standard notation for the nonsingular 
equation is = cu”* with c, m > 0. In this paper the exponent of the nonlinearity is 
written in terms of n = —m with n > 0, and the detailed calculations are done for n > 0. 
The reason for this notation is to avoid the use of negative exponents that might confuse 
the reader in interpreting the result^. The space dimension is N, mostly N = 1 here. 

Theorem 1.1. Let N = 1 and let $(u) = with n > 0 or ^{u) = log(u) (case n = 0). 

Equation ([HD posed in Q = W X {0, oo) with initial data 

(1.3) M(a:, 0) = uo(^) G 7y^(]R), uq > 0 , 

admits a positive very weak solution if s > 1/2 and 0 < n < 2s — 1. There could be 
non-uniqueness of the solutions, but we construct a unique limit solution for every initial 
data, and we prove that it is maximal among all solutions. This solution exists globally in 
time, u G C([0, oo) : L^(M)), and is positive everywhere. 

The range of exponents 1/2 < s < 1 and 0 < n < 2s — 1 is almost sharp. There is indeed 
another isolated case of existence of integrable solutions in ID, namely s = 1/2 and n = 0 
(logarithmic diffusion). The very peculiar properties of this case deserve a separate study, 
but we give a preliminary idea in Section [12] that supports the assertion of existence of 
solutions at least for short times. For the other exponents s € (0,1), n > 0, in dimension 
one, nontrivial solutions do not exist by the mentioned results of [Tj. 

The construction of solutions proceeds by approximation, taking approximate initial data 
that are uniformly positive, so that the problem is no more singular. Passing then to the 
limit in the approximations we obtain a solution that is shown to be non-trivial after some 
effort. It is called the limit solution (upper limit solution, to be precise). It is subsequently 
proved to be a very weak solution. Here we define very weak solution of equation as 

^ By trivial solution we mean u{x, t) = 0 in the whole domain of definition. 

^The reader may also wonder, why the minus sign in the coefficient of It is needed to make an 
increasing function, so that the equation will be parabolic in some sense. 


2 



a nonnegative function u € C{{0,oo) : L^(M^)) such that 


(1.4) 



POO r 

Jo Jr 


4>(u) (—A)®C dxdt , 


and the last integral is absolutely convergent for all C smooth and compactly supported. In 
the theorem 4>(n) = <I>„(u). 

In the course of the present paper we will also establish the main properties of the con¬ 
structed solutions. We select here the main results for easy reference. 

Theorem 1.2. The limit solution preserves mass, f u(x,t)dx = f uo(x)dx. Moreover, it 
is a weak solution for t > 0, and satisfies the bounds 


(1.5) Ci{t){l + |x|2)^/(i+-) < u{x,t) < C 2 \\uo\\it-^ 


with a = l/(2s — n — 1), 6 = 2sa, a continuous function Ci{t) > 0 that may depend on the 
solution, and a constant C 2 {n,s) > 0. The collection of limit solutions generates an ordered, 
-contraction semigroup in L)|_(R) 


A specially important feature of the paper is the construction and properties of the fun¬ 
damental solutions. 

Theorem 1.3. There exists a special function of the form 

( 1 . 6 ) U{x,t)=t-^F{xt-^) 

that is a very weak positive solution of the problem for t > t > 0 and takes on a Dirac 
mass as initial data, u{x,t) —>■ 5(x) as x ^ 0 in the sense of positive Radon measures. The 
profile F is positive everywhere, integrable, symmetric, F{x) = F{—x), and F monotone 
decreasing for x > 0. Moreover, 

(1.7) lim = C{s,n) > 0. 

|ic|^oo 

The constant C{s,n) can be calculated as the constant appearing in the Very Singular So¬ 
lution, a special solution with formula U{x,t) = that has a 

non-integrable singularity at x = 0. Finally, the solution with initial data M 6{x), M > 0, 
is just 

(1.8) Um{x, t) = M U{x, , so that Fm{^) = 


Using the terminology of m we call a fundamental any solution with u{x,0) equal to a 
Dirac delta, and Barenblatt solution a fundamental solution that is also self-similar. On the 
one hand, the Barenblatt solutions of the theorem play an important role in completing the 
existence theory described in Theorem 11.21 In the theory different comparison theorems are 
also used, in particular a new Shifting Comparison result, that we prove as Theorem 14.21 
On the other hand, the Barenblatt solutions explain the asymptotic behaviour of general 
solutions, according to the following general theorem. 
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Theorem 1.4. Let uq € L^(M), let M = f uo(x) dx, and let Um be the self-similar Baren- 
blatt solution with mass M. Then as t ^ oo the solutions u{x,t) and UM{x,t) are increas¬ 
ingly close and we have 

(1.9) lim = 0, 

E —>-00 

Indeed, convergence happens in all norms, 1 < p < oo, in the form 

=0, ^P = p(2s-I+ n) ' 

There are some other results worth recalling. Thus, as a side result of our analysis, we 
construct in Section [ 8 ] the Very Singular Solution (VSS), that is explicit (see Theorem 18. ip 
and has very special properties. Very Singular Solutions have played a special role in the 
theory of fast diffusion equations, as attested e. g. in [ 12 ]. Our VSS will give us a first clue 
to the lower bound for the spatial decay of all positive solutions, that we 

have stated in (lESI) and plays a role in different passages of the existence theory that we 
will develop below. 

After all these theorems are proved and shifting comparison is established, we devote a 
short section to a preliminary presentation of the special case s = 1/2, n = 0. The paper 
concludes with a section on comments, extensions, and open problems. 

Precedents and commentary. (1) Many results are known about Problem ()l.ip - (ll.3p . 
mainly for standard diffusion s = 1, where the Laplacian is used instead of the fractional 
Laplacian. The nonsigular case 4>(u) = u™ with m > 0, is known as the Porous Medium 
Equation when m > 1, the Heat Equation for m = 1 and the Fast Diffusion Equation, 
0 < m < 1 ; their theory has been studied in great detail and is described in monographs 
like O [iTl [Ml [39]. As a basic existence result, each of these equations generate a mild 
solution for every initial data uq £ and the collection of such solutions forms an 

ordered contraction semigroup for every fixed m. 

(2) For equations with fractional Laplacians of the form (jl.ip - (jl.2h with 0 < s < 1 and the 
same of power-like nonlinearity <h(u) = tt™', m > 0, the study of the Cauchy problem with 
nonnegative data in L^(M^), V > 1, has been done in the papers |18l 119] . and most of the 
basic results are still true though the techniques may be quite different. More precisely, 
existence and uniqueness of solutions in the class of very weak or strong solutions have been 
proved, and the main qualitative and quantitative properties are established. Thus, when 
N{m — 1) + 2s > 0 the solutions are positive everywhere in Q = x (0, oo); the so-called 
smoothing effect asserts that initial data produce bounded solutions for positive times 
and indeed 

u{x,t) < C{N, m, s)||no||/ 

where a = N/{N{m — 1) -|- 2s) and 7 = 2s/{N{m — 1) -|- 2s), both positive in this range. 
A main feature of the theory is the existence of fundamental solutions and their use in 
establishing the asymptotic behaviour of general solutions. This was proved in m, also 
under the necessary restriction N{m — 1) -|- 2s > 0. 

(3) Since this condition on the exponents to obtain good behaviour becomes formally m > 
1 — 2s in ID, a similar theory could be expected to hold for very singular exponents m < 0, 
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if 2s — 1 > n = —m > 0 when N = 1. However, the difficulties of dealing the singular 
nonlinearities prevented the inclusion of this extension in the works [TSl [T9] and [TO] . We 
supply in this paper the approach and tools to fill such a gap in the range of exponents 
of Theorem 11.11 In particular, we construct the fundamental solutions and show that they 
are responsible for the asymptotic behavior of general solutions. The text below shows that 
such extended theory is far from immediate and needs some involved tools. 

(4) The case of singular powers ^{u) = —u~^ with n > 0, or <ho(tt) = log(u), was considered 
by the author many years ago in |35j for the standard Laplacian, s = 1. A remarkable result 
of non-existence of integrable solutions was proved for all n > 0 if > 3, for n > 0 if iV = 2, 
and for n > 1 if = 1. That paper is a remote precedent for the present work (previously, 
non-existence for the particular limit case n = 1 in A^ = 1 had been proved in [22] using 
a special transformation). More precisely, when we perform the natural approximation by 
regular problems, as explained in the next section, the sequence of approximations collapses 
to zero for all x G M and all t > 0 for all initial data in the integrable class. This radical 
phenomenon is called instantaneous extinction. On the other hand, existence occurs in the 
remaining cases, A^ = 1, 0 < n < 1, and N = 2, n = 0. This agrees with the results we 
are going to prove in the fractional case. Interesting properties arise in the existence cases, 
see a detailed account in |39l Section 9]. The non-existence results of [35] were extended to 
optimal classes of (non-integrable) initial data in [Ml [151 (IS] • 

(5) In the case of singular nonlinearities and fractional Laplacians, which is our framework 
here, the non-existence of solutions has been recently established in collaboration with 
Bonforte and Segatti, [7], in the range that perfectly complements the positive result of 
Theorem 11.11 plus the announced existence result for s = 1/2, n = 0 in ID. Non-existence 
happens for all singular power cases (n > 0) in dimension N > 2. In all those non¬ 
existence cases we got instantaneous extinction for the initial value problem with any data 
uq G L],_(M). Putting these results together, we obtain a complete picture of the solubility 
problem with integrable data for all singular parameters. 

More on notations. We use the sign / ~ 5 to denote that both functions are proportional 
in a certain limit or range of values I (which may be explicit or understood from the context). 
If the proportionality ratio goes to 1 in some limit then we write f ^ g. In the proofs we will 
often use rearranged functions defined on the line. This means that they are nonnegative, 
symmetric and monotone nonincreasing for x > 0. Other notations will be explained as 
they appear. 


2 Problem, approximation, and limit solutions 

Let us discuss the way to prove existence for the Cauchy Problem (ll.lll - (ll.3l) . i. e., to find 
solutions of the equation posed in Q = M x ( 0 , 00 ) with T > 0 , taking on initial data 
uo{x), assumed to be nonnegative and integrable. In this section we also assume that uq 
is bounded, a restriction that is made for convenience and will be removed later on. The 
nonlinear function is defined, increasing and smooth for tt > 0 , with <h(s) —)■ —00 as 
s —>■ 0. More precisely, we will construct a compete theory for the case where <I> is chosen 
from the list n > 0, mentioned in the Introduction. From this moment on we assume n 
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to be fixed. Following [35j . a strategy of proof of existence or non-existence of solutions is 
based on approximating problem ([LlD-dOl) by the family problems 


dtUe + (-A)^(^>(tt£)) = 0, n > 0, 
tte(x, 0) = uo{x) -I- e for X G 


for any e > 0, so that we avoid data with values on the singular level u = 0. The standard 
theory applies to these problems and a classical solution u^{t,x) exists for all e > 0, and 
it is strictly positive: > s (see details in the next subsection). Moreover, the maximum 

principle holds for these classical solutions and we have > u^i for e > e' > 0. Therefore, 
we can take the monotone limit 


(2.2) tt(x, t) = lim Ue(x, t). 

£—>■0 

This function is a kind of generalized solution of the problem, that belongs to the class of 
limit solutions. It is now an important step of the theory to decide in which sense this limit 
solution is a solution of the equation in a more traditional functional sense (like very weak, 
weak, strong or viscosity solution), and also in which sense it takes the initial data. In 
cases of non-uniqueness of such solutions, the unique limit obtained by the above method 
has been called by various names: maximal solution, SOLA, proper solution,... However, 
this kind of considerations are not the main issue of this paper which is concerned with 
describing the existence and behaviour of the class of limit solutions. In order to recall the 
way the limit is taken (via approximations from above) and to avoid possible confusions, we 
propose the more precise term upper limit solutions for the limits (12.21) . but we will allow 
the simpler name limit solutions when there is no fear of confusion. 

This approximation method has been used in [35] to prove non-existence in the case of 
standard Laplacian as mentioned above, and in [7] to prove the non-existence results for 
fractional diffusion and those singular that do not fall into the cases treated in this paper. 


2.1 Existence and properties of the approximate solntion 


It is convenient to write Us(x,t) = Vs(x,t} -\- e and then try to solve the Cauchy problem 


dtv + (—A)^(<I>£(u)) = 0 with <I'£(u) := $(u + e) — ‘I’(e) 
n(0) = no for x € M. 


for all e > 0. This is a modified problem prepared to avoid the singular level n = 0 of the 
equation by displacement of the axes. Note that for e > 0 and for nonnegative arguments 
$£ is a smooth, positive, monotone increasing function with $£(0) = 0 and ^(.(n) positive, 
bounded and decreasing for all n > 0; ^’^(n) is uniformly positive if v is bounded. The 
theory of existence and uniqueness of weak solutions to the nonsingular Problem (|2.3p 
is given in im Theorem 8.2]. After obtaining these solutions we restore for any e > 0 the 
original u-level by dehning := + e, as stated at the beginning. Clearly, we have that 

Us > e (actually, Us > e) and hence Ug > 0 in Q. 

Let us list some further properties of Vs and Us- For the proof we may again refer to 

[HlllS] and [7]. 
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• Boundedness and regularity. These solutions are shown to be bounded for strictly positive 
times. More precisely, for every t > 0 and every p € [1, oo] there holds 

(2.4) ||u£(-,t)||p < lluollp 


As a consequence, if uq belongs to L°°, then is regular enough to satisfy the equation in 
the classical sense at least when t > 0 , by the results of [l2]. Therefore, = Ve+£ is smooth 
and satisfies the original equation in the classical sense in Q. Under these circumstances, the 
initial data are also taken, at least in the sense of convergence in For unbounded 

data this result about initial data follows from density and contraction in L^, see next 
paragraph. 

• -contraction and comparison. The evolution ()2.3p is an contraction, namely for 
two solutions we have we have 

(2.5) / {uie{x,t) - U 2 £{x,t))+dx < / (uoi - no 2 )+dx for f > 0, 

JRN 


Here, (•)+ denotes the positive part function. In particular, standard comparison follows: 
if itoi < uo 2 a-e., then for every f > 0 we get ui£{-,t) < U 2 £{-,t) a.e. 

• Mass conservation. Nonnegative solutions to the evolution equation (j2.3p conserve the 
mass, cf. [laiig. More precisely, we have for all t > 0 


( 2 . 6 ) 



V£{x, t) dx 



i.e.. 



e) dx 



uo{x) dx. 


• Monotonicity with respect to e. An easy version of the above comparison argument 
shows also that for 0 < e < e' we have 0 < e < Us < up. 

• Time monotonicity. There is an important monotonicity property valid of all nonnegative 
solutions, known as the Benilan-Crandall inequality 

(2.7) dtUe < ^ V(x,t) G Q. 

[n + l)t 

The argument only uses the scaling invariance of the equation and the maximum principle, 
so [5]’s argument applies. 

• The smoothing effect. It says that all solutions with integrable data are in fact bounded 
for positive times. This follows from Theorem 8.2 of [l2], that we adapt to our dimension 
and notations as follows. 

Proposition 2.1. Let $ G C'^(M) be such that <k'(rt) > C\u\ "■ ^ for some n G M and 
|n| > C. If uq G L^(R^) n LP(M^), where p > 1 satisfies 2sp > 1 + n, then a weak L^- 
energy solution to the Cauchy problem (fTTH-lfOl) is bounded m M x (r, oo) for all t > 0 . 
More precisely, it satisfies 

( 2 . 8 ) sup |tt(x,t)| < max{C, Cl 

with Op = l/{2.sp — n — 1) and 5p = 2spyp, the constant Ci depending on n, p, a,C. 
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See also m in that respect. The statement of [42] does not assume that m = —n > 0, 
only that 2sp > 1 — m. The assumptions on <1> are satisfied by the nonlinearities of the 
approximate problems (12.3p . hence the result applies to the approximations Recall that 
we write m = —n and note that the constants in the formula may also depend on e. We 
will solve the latter difficulty later on. 

2.2 Passing to the limit 

We may now pass to the limit e —)• 0 using the monotonicity in e of the family u^. By the 
monotone convergence theorem the limit u is taken in the local sense, i.e., in L^{B) for 
every compact subset B of x [0, oo]. We have 

Proposition 2.2. If uq is a nonnegative funetion in there exists the monotone 

limit u = lime_,.or(e with local convergence in L^{Q). 

3 Existence of nontrivial limit solutions 

The main problem with this procedure concerns the possibility that the limit may become 
identically zero in Qt for some relevant class of initial data. This is what happens for 
V > 2, as established in [7]. In the cases we study here this failure of existence will not 
happen as we will show below. In fact, the limit will be nontrivial for all nontrivial initial 
data. The proof of this general result is long and proceeds in steps. The first of such steps 
consists in exhibiting at least one nontrivial solution. 

A reminder: in the next sections we concentrate on the case of exponents s > 1/2 and 
0 < n < 2s — 1. Note that even if the diffusion is singular, the range is formally the same 
as the good fast diffusion range 1 > m > {N — 2s)/N, considered in the general theory of 
m, hence it is supercritical in the notation of that paper. But there only exponents m > 0 
were considered. 


3.1 The very singular solution, I 

In paper [40] a formal solution of equation ut + (—= 0 is constructed with the form 

(3.1) U{x, t) = H{t)F{x) = C{N, s, m) \^\- 2 s/(i-m) 

when m > [N — 2s)/V, so that the spatial profile has a non-integrable singularity at x = 0. 
This type is called very singular solution in the literature (VSS for short). It is also proved 
that such formal solution is a limit of a monotone increasing sequence of standard solutions. 

Following that paper we try the same formula here for m = — n < 0 with factors 

Hi(t) = Fi(x) = 

• Let us check that it works in the parameter range, s > 1/2, n + 1 < 2s with n > 0. We 
have 4>{F{x)) = — so that a simple calculation gives 

(3.2) C^4>{F{x)) = -K{s, n) |a;| 2 W(i+n)- 2 . ^ _ 


We prove in the note below that K[s, n) is positive in this range of values of s, n. It is now 
easy to see that if we put 

(3.3) Ui{x,t) = Hi{t)Fi{x) with = K{n + 1), 


then C/i is a solution of the ID equation in the range 0 < n < 2s — 1, unless at the singular 
point X = 0 . 

The singularity prevents the VSS from being acceptable as an example of nontrivial so¬ 
lution in the theory we are considering for Problem (ll.ip - (jl.3l) . at this stage. But we will 
return to this topic in style in Section 0 

Note on fractional Laplacians of power functions. The s-Laplacian of a power, 
(-A)I x|", a > 0 , a 7 ^ 2 s, is the power k{a,s) |x|“ with a constant factor that for 
= 1 equals 


(3.4) 


k{a, s) 


, 25 r((l + a)/2) r((-Q; -I- 2s)/2) 
r((l + a- 2 s)/ 2 )r(-a/ 2 ) 


In our case we take a = 2sn/(l -|- n) > 0. Hence, r((l -|- a)/2) > 0; besides, (—a + 2s)/2 = 
s/(l -|- n) >0 so that r((— a -|- 2 s)/ 2 ) > 0; moreover, a/2 = sn/{l -|- n) < 1 (since s < 2 
and n/(l -|- n) < 1/2), so that r(— a/2) < 0. Finally, for our choice of a 


1 a — 2 s 


2ns 2sn — (1 -|- n)( 2 s — 1 ) 

- + l-2s = -^- - - - 

1 + n 1 -|- n 


1 -|- n — 2 s 
1 -|- n 


G(0,1). 


It means that r((—a -|- 2s)/2) < 0. Therefore, k{a, s) > 0 for theses particular values of a 
and s. This is what we have called K{s,n) some lines above. 


3.2 A bounded subsolution 

Though the singularity prevents the VSS from being directly useful as an example of non¬ 
trivial solution, the idea is not completely lost. We will start from it to construct a useful 
smooth variant, but it will be only a subsolution. 

Lemma 3.1. There exists a bounded and smooth function in separate-variables form 
(3.5) U{x,t) = H{t)F 2 {x) 

which is a positive subsolution of the equation for some 0 < t < T and all x G M. Moreover, 
F 2 is symmetric, radially decreasing and F 2 {x) ^ as x ^ 00 . Changing the 

form of H we can get a supersolution. In both cases H is continuous and the initial value 
H{0) >0 can be chosen. 

Proof, (i) In order to avoid the problem with the singularity of the VSS, we round the 
function Fi in a small ball near x = 0 to get a smooth positive F 2 , so that instead of the 
exact formula {—Ay(p{Fi) = —K{s)Fi{x) for x 7 ^ 0 we get an approximate equation for 
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all X that we can still use. Indeed, in view of the perturbation we get that {—A.y(f>{F2) is 
bounded on bounded sets; on the other hand, the difference 

(-A)"(/.(Fi) - {-Ay<P{F 2 ) = 0(|xr(^+2'*) as |x| ^ oo. 

The last formula comes directly from the representation formula for (—A)^ plus the fact 
that 4>{Fi) — 4>{F2) has compact support. This correction 0(1x1“^^"*“^®^) is lower order with 
respect to {—Ay4>{Fi) at infinity since (—A)'^(^(Ti) ~ Ti = T 2 ~ for all large 

|x|. Summing up, there exist constants Ki,K 2 > 0 such that 

(3.6) - KiF2ix) < -(-A)XF2) < K2F2{x) . 

If we now take H{0) = a > 0 and H' < we get a subsolution, 

dtU + {-Ay^n{u)<o, 

in some time interval. In the last case maybe H decays and vanishes in finite time (even 
if this is not realistic for actual solutions, as we will see, it is just the form of the modified 
subsolution). Take F = F 2 and this choice of H to end the construction of the desired 
subsolution U. Finally, note that since F[{0) can be taken at will we find a family of 
bounded subsolutions and the norm can be taken as small as wanted. 


(ii) In the same way, if we take Fi as before and F[(0) = a > 0 such that F[' > K 2 H ”, 
then we get a formal supersolution. In this case it will exist for all times. □ 

• Case n = 0. This case is settled by replacing the power —by log(u) and repeating 
the above procedure. We will need the calculation the s-Laplacian of the logarithm. We 
have 


(3.7) 


(—A)®(log |x|) = c(s) |x| 


A short proof is as follows: using the previous formula for a > 0, a very small we get 


(3.8) 


k{a, s) 


,2,_2 (2g-i)r(i/2)r(s) 
r((3 - 2s)/2) 


where we have used r(— a/2)) ~ r(l)(— 2/a), and r((l+a— 2s)/2) ~ 2/(1—2s) r((3— 2s)/2). 
We now use the expression log(x) = lima_j.o(x“ — l)/a, x > 0, to conclude that the formula 
is true with c(s) = limc^o k{a, s)/a > 0 if 2s > 1. 

Then the rest of the steps is quite similar and the conclusions of Lemma l3.II hold. □ 


Proposition 3.2. Let N = 1 and 0 < n < 2s — 1. Let the initial data uq he positive 
and integrable and satisfy uo{x) > CF 2 {x) for all x. Then the limit solution constructed 
as in Section [2] is non-trivial. In fact, it sits on top of one of the constructed subsolutions, 
hence it is a positive very weak solution of the equation, at least in a certain time interval 
0 < t <T. Nontrivial solutions are also obtained for n = 0, 1/2 < s < 1. 
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Proof. The subsolution can be compared by classical results with all approximate problems 
and remains below for all the existence time. This allows us to prove that the limit solution 
does not vanish identically when we take initial data that decay equal or slower than F. 
By the contraction the same nontrivial limit happens for any continuous and bounded 
initial data. 

Remark. This is the first successful step in a long road that leads to the proof that all 
upper limits corresponding to nontrivial data are in fact positive weak solutions. 

3.3 Some properties of nontrivial limit solutions 
• Scaling property. If u{x,t) is a nontrivial limit solution, then so is 
(3.9) UAL{x,t) = Au{Lx,L‘^^A~^^^"'h) 

for all parameters A,B > 0. We leave the easy proof to the reader (see similar arguments 
in [39]). 


• The smoothing effect. It says that all limit solutions with integrable data are in fact 
bounded so that we use bounded solutions in the proofs. We have for all limit solutions 

(3.10) sup|n(x,f)| < C' 2 f“"||uo||f 


with a = l/( 2 s — n — 1 ) and b = 2 s/( 2 s — n — 1 ), the constant C 2 depending on n, s. 

Proof Use the smoothing effect in the rough form already proved for the approximate 
solutions and let t = 1 to conclude that the result holds for M = ||rto||i = 1- When M 7 ^ 1 
and t ^ 1 use the scaling rule (j3.9p in the usual way to reduce the proof to the particular 
case. This is a well-known scaling trick. □ 

• Time monotonicity. The limit solutions satisfy 


(3.11) 


dtu < 


u 

(1 + n)t ' 


It follows from the same property for the approximations. As we said there, the argument 
has a proof using scaling arguments originally due to Bmilan and Crandall [5]. 

Also the properties of boundedness, pointwise comparison, and contraction of the 
approximate solutions pass to the limit without change. 

• Space monotonicity. Aleksandrov’s principle. The Aleksandrov-Serrin reflection 
method is a well-established tool to prove monotonicity of solutions of wide classes of (pos¬ 
sibly nonlinear) elliptic and parabolic equations, cf. [Ill32|. It has been quite useful in 
particular in the case of the PME, as documented in [IIIEHI. This is the version proved in 
m Theorem 15.2] for nonlinear parabolic equations with fractional diffusion of the type 
dm), and adapted to our situation 

Proposition 3.3. Let the unique solution of (12.3|] with initial data uq € L^(M), uq > 0. 
Under the assumption that 


(3.12) 


1 ^ 0 ( 3 ;) < uo{2a — x)) for x > a 
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for some a > 0, we have for all t > 0 

(3.13) V£{x,t) < Vs{2a — x,t) for x>a. 

In plain words, the result deals with comparison of a solution with its space reflection with 
respect to the point x = a. If it is true for t = 0, then it is true forever. After passage 
to the limit the same comparison is true for the obtained upper solutions. An immediate 
consequence of this that has been used in the literature and we will use below is 

Corollary 3.4. An upper limit solution with initial data supported in the half-fine {x < 
a} is monotone nonincreasing in x in the region {x > a, t > 0}. If the initial data are 
supported in the half-fine {x > —a} the solution is monotone nondecreasing in x in the 
region {x < —a,t > 0}. 


4 Comparison results 

The standard comparison theorem (Maximum Principle) applies the approximate prob¬ 
lems, hence it will be valid in the limit e —0 for the class of upper limit solutions. The 
Aleksandrov principle is another comparison theorem. In the sequel we will use two other 
comparison results, that we discuss next. 

4.1 Symmetrization and concentration comparison 

Symmetrization techniques are a very popular tool of obtaining a priori estimates for the 
solutions of different partial differential equations, notably those of elliptic and parabolic 
type. The application of Schwarz symmetrization allows to obtain sharp a priori estimates 
for elliptic problem by comparison with a model symmetric problem. For parabolic problems 
the usual pointwise comparison of the solutions of the two problems fails, and is replaced 
by comparison of integrals, [3]. In the case of the porous medium equation ut = Au^ that 
result was established in I331E7], and holds for all m > 0. In order to state the result we 
will use, the following definition is needed: 

Definition. Let f,g ^ be two radially symmetric functions on . We say that 

/ is less concentrated than g, and we write f ^ g, if for all i? > 0 we get 

(4.1) / f{x)dx< / g{x)dx. 

iBjiiO) JBjiiO) 

The partial order relationship ^ is called comparison of mass concentrations. In the 
applications we are going to assume that / and g are rearranged functions. 

The following result is proved in |43j . 

Theorem 4.1. Let Let ui,U 2 be two nonnegative, weak solutions of the equation ut -f- 
(—A)^<I>(u) = 0, posed inQ = x (0, oo), with nonnegative initial data uqi, uo 2 G L^(M^). 
Assume that both uo 2 and uoi are rearranged and uo 2 ^ uqi. Assume moreover that the 
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nonlinear function $(w) is positive, smooth and concave for u > 0. Then, for all t > 0 the 
functions ui{-,t) and ui{-,t) are rearranged and we have 

(4.2) U2{-,t) ^ui{-,t). 

In particular, we have ||tt 2 (-, i)||p < ||ui(-,t)||p for every t > 0 and every p G [l,oo]. 

This result applies to the solutions of the regularized problems (12.11) . In the limit it 
will apply to all nontrivial solutions of the singular fractional FDE (11.11) with and 

n > 0, for all 0 < s < 1. 

Note that in ID the integrals in formula (j4.1l) have a more classical interpretation: the 
integral f{x)dx is just the distribution function of the mass density /, that we are 
assuming to be nonnegative, integrable and monotone decreasing for x > 0. Therefore, 
formula (14.2p is just a comparison of distribution functions. 


4.2 Shifting comparison 

A new comparison result that is related to symmetrization in spirit and techniques is based 
on lateral displacement of the solution, viewed as a mass distribution. 

Theorem 4.2. Let two functions uoi )^02 S T^(M) with the following properties: 

(i) They are nonnegative and rearranged around their points of maximum xi and X 2 resp., 
with xi < X 2 - The total mass is the same. 

(a) We assume moreover that uo 2 dx < uqi dx for every x G M. 

Besides, we assume that is monotone, concave and defined on R+ with 4>(0) = 0 and 
0 < $'(0) < oo. Then, the following comparison inequalities hold for the corresponding 
(limit) solutions 

/ X 

U 2 {x, t) dx < 

-oo 

for every x G R and every t > 0. 

This result is called the Shifting comparison lemma. It is essentially one-dimensional and 
it was established in the PME case by Vazquez [M] and it proved useful in studies of free 
boundary location or asymptotic behaviour. It is related to mass transport and Wasserstein 
distances, |45l 136] , It will be crucial in some proofs below, like the proof of the existence 
for general initial data and the asymptotic behaviour. Since it has a rather technical and 
long proof, we will delay to Section [TT] at the end of the paper. 



5 Mass conservation and global solutions 

The property of mass conservation plays an important role in passing from local-in-time 
existence to global solutions, since it prevents the phenomenon of finite-time extinction. A 
first result is as follows. 
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Proposition 5.1. Let us assume that N = 1, s > 1/2 and 0 < n < 2s — 1 and let us 
assume that uq is integrable and uo{x) > cF 2 {x), i.e., 

(5.1) uo(a:) > c/(l + 


for some c > 0. Then the limit solution not only is positive in a certain time interval 
0 < t < T, but it also conserves mass in that interval: 


(5.2) 




dx. 


Proof, (i) Let us first assume that n > 0. We will prove conservation of mass for small times. 
We take a nonnegative non-increasing cut-off function (/(s) such that C('S) = 1 for 0 < |s| < 1, 
C(s) = 0 for |s| > 2, and define Cr{x) = C(|x|/i?). We have (—A)^(/i € n 

moreover, |(—A)^Ci| ~ as |x| —>■ oo. The radial cut-off function (^r has the scaling 

property 

(5.3) (-A)\r(x) = R-^^i-AYCiix/R). 


We take the approximate solutions of equation O) with nonlinearity as in Section 
[21 multiply by Cr and integrate by parts. We have 


{Ue{ti) - UQ- e) Cr dxdt = 


"■ {-AyCRdx. 


Note that the last integral is absolutely integrable. Passing to the limit e —>■ 0 we get for 
every t > 0, 

(5.4) [ u{ti)CRdx- f uoCRdx = f [ u~'^ {-A)XRdxdt. 

Jr Jr Jo Jr 

Let us split the right-hand side of (15.4p into the integrals for |x| < R and |x| > R. We only 
estimate the integrals in x, forgetting for the moment the time integration. We get 


1^1 (^) I < 


|a;|>R 




|a;|>_R 


Putting X = Ry we get 


h{R) < 


(Jg2ns/{l+n) 



dy 


CR-"! [ 
J\y\>i 


dy 

|y|2+7 


Since 7 = 2s — 1 — 2ns/(l -|- n) = 2s/(n -|- 1) — 1 > 0 , we get Ii{R) —>■ 0 as i? —>■ 00. On the 
other hand, for the analogous integral in the set |a:| < R we get 

\I 2 {R)\ < [ u-^{t) |(-A)^Cr| dx<^ [ |x|2W(l+n) dx < CR-"' 

J\x\<R P J\x\<R 


that goes also to zero as i? ^ 00. Since these estimates do not depend on t (for small t) we 
may go back to equation (15.4h and let i? —>■ 00 to get 




dx, 
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which is the mass conservation law. This law holds for the small times for which we have the 
lower estimate for the limit solution u (for instance, the estimate coming from comparison 
with the subsolution constructed in Section [3l as used in Proposition I3.2p . 

This method of proof is inspired in a common technique that was used to prove the property 
in the case s = 1 for m > {N — 2)/N, N > 2. The beginning is the same, but the further 
details are quite different. 

(ii) In the case n = 0 the nonlinearity is logarithmic and not power-like. We arrive at 
formulas Ii{R) and l 2 {R) with the bound clog(l + \x\^) instead of c\x\^'^^ /and we put 
n = 0 in the rest of the places. The remaining steps follow easily. □ 

Proposition 5.2. If the initial data uq is a rearranged function with mass M > 0, the limit 
solution u{x, t) exists globally in time, is positive everywhere and the mass is conserved for 
all times. 

Proof, (i) We first prove that for rearranged initial data in the same class uq > cF 2 {x) the 
solution exists for all times and conservation of mass holds also for all times. We know that 
u is positive and conserves mass for 0 < f < Ti = Ti{u). Assume that this maximal time is 
finite. We use the scaling property of the equation to define of new solutions 

(5.5) UL{x,t) = Lu{Lx, 

with L > 1. Let us introduce the notation ul = Tlu for future reference. This scaling, 
a particular case of (13.9p . keeps the mass of the solutions u and ul identical. Now, since 
the time of existence of the family ul shrinks to the time Tl = Ti /we seem 
have a problem in using rescaling. But the problem can be fixed by using symmetrization 
(concentration comparison) and we end up with an expected gain. 

Indeed, the new solution for L > 1 is clearly more concentrated than ui{x,t) = u{x,t) 
at at time t = 0, hence it will be more concentrated at all times by the result of Subsection 
O Then, the concentration relation (14.2p immediately implies that the mass of ul must 
be conserved as long as the mass of u is, say until Ti{u)] the justification for Tl < t < Ti 
is done by doing symmetrization comparison on the approximate problems and passing to 
the limit, that cannot be trivial because of this argument. This means that the nontrivial 
existence time with conservation of mass for ul is T{ul) > Ti. Undoing the scaling we get 
the same property for n in a time T{u) > hence Ti must be infinite. 

(ii) Next, we prove that the solution is positive everywhere. Let us do an analysis of what 
happens if a rearranged solution touches zero, and show that this cannot happen. First, we 
use the monotonicity properties in space and time to show that, if the solution vanishes at 
t = ti and X = R, then we must have u{x, t) = 0 for all |x| > R and t > ti. Using the proof 
of the mass formula of Proposition 15.11 we have for the approximations Ue 

4: [ Ueix,t)C{x)dx = - f ^n{ue){-i^yc,dx 

ul Jr Jr 

Taking a cutoff function supported in [—with Ri < R, we know that —(—A)®C > 0 
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for |x| > Ri, and since ^ni^e) tends to minus infinity in the set we have described, we have 


- ^nius) {-Aycdx ^-oo 

J |a;|>i?i 

as e —>■ 0. The integral for \x\ < Ri has a bounded limit since the limit u is bounded below, 
hence |<h(u)| is bounded. Applying this between times ti and ti + r, we can show that the 
weighted mass must be zero for all times larger than ti. Since the solution is rearranged, 
this implies that the mass will be zero, which is excluded by the previous step. □ 

(hi) The next step is to consider rearranged initial data that do not satisfy a bound from 
below like (EtD, for instance uq may be compactly supported. Let M > 0 be the initial 
mass. In that case we make a small perturbation by adding to uq a tail of the form 
r]^{x) = h(l + and besides we truncated the initial data on top to make it 

bounded. We may do all this and conserve the mass M. In this way we obtain a solution 
to the problem where the previous analysis applies, and mass is conserved. We now use 
the contraction property and the conservation of mass for u^, we conclude that the limit 
solution u corresponding to initial data mq niust have mass 



dx > M 


11^0 “ 


which is positive for 8 small, hence tt is a global solution. The argument is justified as limit 
solution, i. e., in the limit of approximate problems. Finally, by letting <5 —)■ 0 we derive the 
mass conservation property for u. □ 


Positivity can now be obtained for all solutions with continuous initial data, not necessarily 
rearranged, by standard comparison with a solution with compactly supported rearranged 
data (after possibly a space displacement to fit it under uq). Such a subsolution is positive 
for all positive times, hence u is too. We will return to this question later on, after we get 
some quantitative estimates on the behaviour. 


5.1 Concept of solution. The mass function 

As a preliminary for the next developments, we need to clarify the type of solution that we 
get at this stage when we pass to the “limit solutions”. In the end, we would like to prove 
that our positive solutions are very weak in stated sense that 


(5.6) 


u Ci dxdt 


^n{u) {—Ayedxdt 


for all smooth test functions with compact support. However, the last term offers a difficulty 
as long as we do not know the decay of u at infinity, i. e., as long as we do not control the 
growth of since typically (—A)®C behaves like 0(|x|^'''^^) as |x| —oo. 

The conditions of Proposition 15.11 do allow for a correct passage to the limit e —>■ 0 in 
the definition, but the conditions of Proposition 15.21 do not. Keeping the assumption of 
rearranged data, we find a remedy by weakening the definition by integration in space to 
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get a new function V (r, t) = u(x, t) dx so we have Vy = u for all r > 0 and we can write 
the equation formally as 

(5.7) Vt= f utdx = - 

Jo 

Putting (—A)^ = —d'^^{—Ay' with s' = 1 — s > 0, and integrating we get 

(5.8) Pi = ap-A)-^'(^>„(n)), Vr = u. 

This is an integrated version of the weak solution that makes perfect sense for the ap¬ 
proximate problems and for their limits in the very weak sense. We will call V{x,t) the 
mass function; it is the distribution function in Probability, but that name might lead to 
confusion here. 



6 The Barenblatt solutions in ID 

After establishing mass conservation for all times we can construct the Barenblatt solutions 
and derive the main properties. 

6.1 Existence 

Take one of the rearranged solutions ui{x,t) of the previous section (Proposition [52]) with 
initial data rtoi(a^) that we may assume continuous. It exists globally in time, is positive 
everywhere and conserves mass. Let us fix the norm of ui to 1. Take then the rescaled 
family {uL{x,t) = Tlu : L > 1} defined by formula (15.5p of the previous section. Finally, 
pass to the limit 

(6.1) U{x,t)= lim UL{x,t). 

L—>-oo 

We have to show that this limit exists and has the desired properties. 

(i) In principle, the family ul{x, t) converges weakly in for t > t > 0, and there may also 
be a non-unique limit. It is best to use an argument based on concentrations, that makes 
it natural to argue with the family of mass functions VL{x,t) = Jq UL{x,t) dx. By the 
concentration comparison result we see that the family {V^} is monotone increasing in the 
parameter L for x > 0 (resp. negative and decreasing for x < 0). Hence, the limit Voo{x,t) 
exists and does not depend on subsequences Lk —oo. The convergence VL{x,t) Voo{x,t) 
is uniform convergence in x for every fixed positive time. 

(ii) Differentiation in x gives the unique weak limit U = limx,_>.oo Ul- Stronger convergence 
will be proved later on. 

(iii) For t = 0 we have ^00(2;, 0) = 1/2 for all x > 0, ^00(2;, 0) = —1/2 for all x < 0, in other 
words the limit of the initial data is a delta function. We have to show that for all positive 
times U{-,t) is not trivial away from x = 0 (i.e., it is not equal to 5{x)). This follows from 
the smoothing effect (proved in Subsection 13.3p that applies uniformly to all functions ul- 
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We conclude that U{-,t) is a bounded function for all t > 0, therefore is not a Dirac 

delta, and Voo{-,t) is Lipschitz continuous in x uniformly if t > r > 0. 

(iv) Let us prove next that the limit must be a self-similar function. The argument is 

based on passing to the limit in the scaling family using the group of transformations. In 
fact, for all L, fe > 0 we have TLTfcUi = TkiUi, Passing to the limit fe —>■ oo and using the 
the uniqueness of the limit because of the comparison of concentrations, we get U{x,t) = 
TLU{x,t) = LU{Lx, , hence LC/(Lx, = U{x,l). In the usual way it 

follows that 

(6.2) U{x,t) = t~°‘F{xt~^), a = l/(2s — (1 + n), 
with F{x) = U{x, 1). 

(v) It is immediate that the profile F is positive everywhere, bounded, integrable, and 
rearranged (i. e., F{x) = F{—x) and F monotone decreasing for x > 0). 

(vi) To construct the Barenblatt solutions with initial data M (5(x) with any mass M > 0 we 

use another scaling T' defined by {T'u){x, t) = Mu{x, that transforms solutions 

of mass 1 into solutions of mass M. We get a self-similar solution with the same formula 
as before, but now the profile is 

(6.3) Fm(x) = . 

6.2 Alternative approach and better convergence 

The idea is to apply the same type of approach to the approximate problems (12.ip to obtain 
a fundamental solution of each of those problems. After careful inspection, we see that we 
can find a solution Uooe{x, t) with initial data C/oe(x, 0) = 6o{x) + e. This is better done by 
using the formulation Ve{x, t) = u^{x, t) — s that solves Problem (12.31) . to which we can apply 
the usual theory and comparison of concentrations. Therefore, we obtain a fundamental 
solution with this argument, the difference is that this time we cannot conclude that it 
is self-similar. On the other hand, the whole collection of rescaled solutions UL,e{x-,t) are 
uniformly bounded for t > r > 0 hence, by the regularity results of [l2], they are uniformly 
(7“ continuous for a certain a > 0. The family is thus locally compact in (in both 
space and time), which means that the convergence UL,e{x-,t) —>■ Uoo,e{x,t) takes place in 
the strong sense of L^(M) for every t > 0. This is a key improvement in the situation. 

Now we take the monotone limit of these fundamental solutions to get 

lim Uoo,eix, t) = Uooix, t) 

£->■0 

By a simple comparison, Uoo{x,t) > U{x,t), where U is the previously constructed Baren¬ 
blatt solution, formula (16.2p . By the equality of masses we conclude that both are the same 
function, Uoo = U. 

It is now easy to see that UL{x,t) —>■ U{x,t) in L^(M) for every t > t > 0. In fact, 

\\{UL{x,t) - U{x,t))+\\Ll(-R,R) < \\iUL,e{x,t) - C/(x, t))+ 

< \\{UL,e{x,t) - Uoc,e{x,t))+\\Ll(-R,R) + o{e), 
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where the last term is the contribution of \\Uoo^e{x,t) — Fixing a small 

e > 0 and and using the convergence UL^e{x,t) —>■ Uoo,eix,t) in L^(M), there is an > 
0 such that have ||C/£,(x,t) — [/(x,t))|| < e for L > L^. All together, this means that 
{Uiix^t) — [/(x,t))+ —7> 0 in L^{—R,R) as L ^ oo. Since the mass in the far field is 
uniformly small, and there is total mass equality, we get — [/(x,t) —> 0 in 

This is the improved convergence that we needed. 

6.3 Uniqueness 

Proposition 6.1. The fundamental solution is independent of the rearranged function uq 
that starts the construction. 

Proof. For bounded functions of compact support this is done by symmetrization compar¬ 
ison and iterated limits. Indeed, let us start from functions ui and ui and let U and U the 
corresponding Barenblatt solutions obtained in the limit of the scaling process. Then it is 
not difficult to see that the scaled function uok is more concentrated than uq if k is large 
enough. Therefore, UkL is more concentrated than uol, hence in the limit U{-,t) is more 
concentrated than U{-,t). The reverse relation also holds, hence U = U. 

For general data we use contraction. □ 


7 Study of the profile 

We want to know more about the profile F of the Barenblatt solution, in particular its 
equation and its asymptotics for large |x|. 

7.1 Profile Equation 

We can apply to U the concept of solution of last section and perform the computation 
when U is selfsimilar. Then, with ^ = xt~°‘, 

Vt = -^[ / {F{^)+CF'{^))dx = -at-^^F{^). 

dt Jo Jo 

Therefore, equation (|5.8I) becomes at t = 1 

(7.1) {F~'^{x)) = ax F{x), x > 0. 

This is the integro-differential equation satisfied by the profile F. Note that 1—2s' = 2s—1 > 
0, so that the operator in the LHS has a positive degree of differentiation, dxC~^ = Fl£J~^l’^. 
We can also write the equation as 

nr 

(7.2) / {-AyF-^{r)dr = axF{x), 

Jo 


F (r) dr = a{—A) ^{xF{x)) 
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7.2 Asymptotic behaviour of F 


According to the analogy with the Barenblatt solutions constructed for the case m > 0 in 
m, since we are in the case m < mi = N/{N + 2s) (a value that was important in that 
respect), we expect the following behaviour. 

Proposition 7.1. The profile F decays as \x\ —>■ oo like r~^, with 7 = 2s/(l + n) > 1. 
Indeed, there exists the finite positive limit 

(7.3) lim F{r) r'^ = , 

r —>-+00 

and moreover F{r) < CooT~^ for all r > 0. 

Remark. Later on we will calculate the constant Coo explicitly as a corollary of our work 
on very singular solutions, Section [ 8 j 

Proof, (i) We need a first estimate of the decay of F. Since it is an integrable and rearranged 
function we immediately get F{x) < C/lxl, which is too rough. A better estimate is obtained 
as follows: We start the scaling procedure to construct U by taking as initial data uoi the 
function F 2 {x) of Lemma (3.II By comparison with the subsolution of the Lemma we have a 
decay rate ui{x) > for all |x| > 1 uniformly in some time interval 0 < t < Ti. 

We can show a similar decay rate for the limit solution U{x,t) by applying concentration 
comparison. Indeed, the mass of ui in region |x| > i? > 1 is estimated as 
and the mass of U has to be less than that by the comparison. By using the monotonicity 
of U w.e.t. \x\ we get the conclusion that U is less than ci(t) . Recall that F{x) 

is just U{x,l). Therefore, 

F{x) < . 

The reader can find similar arguments in @01 Section 12]. 

(ii) The power-like bound from below can be obtained following the ideas of @0]: We 
start from the homogeneity estimate (13.111) that says that {1 + n)tut < u. In terms of the 
self-similar profile, this just means that —(1 -|- n)ot{F + rF'{r)) < T, r = |x| > 0, hence 

—rF'(r) 1 2s 

- — < 1-1 -= -. 

F{r) (1 -|- n)a 1 -|- n 

Integration of this inequality gives the following lower bound, valid for all r > 1, all s € (0,1): 

(7.4) F(r) > . 

Moreover, the function J(r) := T(r) is monotone non-decreasing with r, so that 

it has a limit as r ^ 00 . Let us put 

(7.5) lim F(x) = Coo(s, n). 

|fr|—>-oo 

In principle the limit may be finite and positive or infinite. By part (i) it is finite. 
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8 The very singular solution 

The a priori bounds on the profile of the Barenblatt solutions make it easy to pass to the 
limit M —>■ oo and obtain a special function, called the very singular solution (VSS), much 
as we have done in [30] for m > 0. This is the result 

Theorem 8.1. The Very Singular Solution calculated in Section\^ 

(8.1) U{x, t) = C(n, s) 

is the limit of the Barenblatt solutions UMix,t) as the mass M goes to infinity. 

Proof. We recall that the value of the constant C(n, s) can be explicitly computed from the 
calculations there as C{n, s) = K{n, s)(l + n). 

By the established comparison properties, it is clear that the sequence of Barenblatt solu¬ 
tions Um{x, t) is monotone increasing with M > 0. Next, we check that they are all bounded 
above by the VSS. A direct comparison of Um with U is difficult to justify directly, hence 
we argue in another way to get our conclusion. 

We first fix the mass equal to one and use the upper profile bound, F{r) < 
for all r > 0, to conclude that 

Ul{x,t) < . 

We know that the whole sequence Um can be obtained from M = 1 by the rescaling 
UM{x,t) = {TMu)ix,t) ■= M u{x, We immediately see that Uuixfi) satisfies 

the same upper bound, even with the same constant, Coo(M) = Coo(l)- Once, we have the 
same upper bound for the whole sequence, we may pass to the monotone limit and get 

Voo(x,t) 

But since the functions Um are invariant under the mass conserving scaling, so is the unique 
limit, hence Uoo{x,t) is self-similar, Uooix,t) = t~^Foo{xt~°‘). Also Foo{y) > FM{y) for all 
y > 0 and all M > 0, and Foo{y) < It easily follows from a tail analysis 

that Fooiy) = hence 

TT (.j. „ ^l/(l-|-n)| |-2s/(l-|-n) 

^ 00\'^ ^ ^ J — '~'00 ^ 1*^ I 

We conclude that Uoc is another possible very singular solution of the equation obtained as 
limit of upper limit solutions. In order to see that Uoo and U must be the same we only 
have to check that both are weak solutions for the equation for x 7^ 0, in other words, the 
profiles must be weak solutions for the F equation dZI]). This is what selects the constant 
in a unique way. We conclude that Coo = C{n, s) and Uoo = U. 

As a corollary of the above result we can refine the information on the Barenblatt solutions 
obtained in Proposition 17.11 as follows 

Proposition 8.2. The profiles Fm, M > 0, of the Barenblatt solution satisfy the uniform 
bound F{r) < C{n, for all r > 0 and this estimate is sharp at infinity 

(8.2) lim F{r) = C{n, s), 

r—^+oo 

where C{n,s) is the constant of the VSS. 
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9 The theory for general initial data 


Here we want to complete the proof of Theorems II.11 and 11.21 on the existence and properties 
of solutions of the Cauchy Problem. We establish existence, positivity and behaviour as 
|x| —>■ oo. The basic analysis is done for compactly supported and rearranged data. Then 
we perform a series of extensions of the results to greater generality. 

• Compactly supported and rearranged data, (i) Existence and conservation of mass. 
Assume to fix ideas that uq is supported in the interval [—R, i?], and has total mass M. We 
may use shifting comparison with the Barenblatt solutions with the initial masses located 
on either end of the support to prove bounds from above and below for the mass function 
of u, defined by v{x, t) = f_^ u(x, t) dx, in terms of displaced versions of the mass function 
of the Barenblatt solution with the same mass, Vm corresponding to Um- We get 

Vm{x — R,t) < v{x,t) < Vm{x + R,t). 

Recall that VM{oo,t) = M for every t > 0 by mass conservation. This not only proves 
non-trivial existence for all times, but also conservation of the total mass and existence of 
a nontrivial mass (i. e., integral) on every interval of length larger than 2R. 

We have to justify that shifting comparison applies to Barenblatt solutions, and this is 
done by starting the construction with compactly supported ui{x) and its scalings uq^l{x) = 
Lui{Lx). 

(ii) Positivity. By the Aleksandrov principle, see Corollary 13.41 the function u{x, t) is 
monotone decreasing in x for x > R and t > 0. Together with the previous mass analysis 
this implies that u{x,t) > 0 for x > R. A similar argument happens for x < —R and 
ensures for positivity for all \x\ > R and all t > 0. 

In order to establish the positivity for |x| < R at times t > 0 we take time ti > 0 and 
move the the origin or coordinates to a point xi > 2R. Setting y = x — xi we see that 
u{y + xi,t) is positive near the new origin and bounded below by a rearranged function 
uo{y) with small support [—to which the preceding result applies so that u{y,t) is 
positive for \y\ > Ri. By comparison u{x,t) is positive with a uniform lower bound for 
|x| < R and t = 2ti. 

(hi) Asymptotic behaviour as |x| —>■ oo. 

Proposition 9.1. For every solution with rearranged and compactly supported data with 
mass 1, we have 

(9.1) 

for all \x\ > 2R, and the constants do not depend on the particular u, nor on t. 

Proof The upper bound comes from the mass analysis of (i) and the monotonicity in x for 
large |x|. We get 

nx POO 

X u{x, t) <2 I u{x, t)dx <2 / tt(x, t) dx 
Jxl2 Jxl2 
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and by the Shifting comparison resnlt, this mass is less than the mass of the displaced 
Barenblatt Uoo{x — R,t) which is proportional to 

7 = (2s - 1 - n)/(l + n). 

We conclude that there is a constant ci such that for all x > 2R. 

• The lower bound follows from similar arguments, but now we compare with the mass of 
Uoo{x + R,t). We get 

rkx roo roo 

{k — l)xu{xR)> / u{x,t)dx= / u{x,t)dx— / u{x,t) dx = R — I 2 ■ 

J X J X J kx 

The first integral Ii is estimated from below by the mass of Uoo{x + R,t) in the same 
interval which is accurately given by , while the second is estimated from 

above by the mass Uoo{x — R,t) and gives \kx\~'^ in first approximation. Hence, 

for k > C'2/Ci and x large enough we get u{x,t) > which implies the 

stated lower bound. 

Remark. The translation of these results for data with mass M 7^ 1 is easy by using the 
mass-changing transformation T. 

• Rearranged data. If the initial data are rearranged but not compactly supported, 
we use approximation of the data from below with compactly supported data, so that by 
comparison the property of positivity follows. Conservation of mass comes from the 
contraction property. The asymptotic lower bound in (19.111 still holds, but the upper bound 
need not hold (it depends on the behaviour of the initial data for large |x|). 

• Continuous data. In this cases there is a possible problem with the positivity of the 
solution, that is still not guaranteed by us. This is proved for continuous data uq by putting 
below some rearranged data with compact support, uq, to which the previous theory applies. 
Then we can apply comparison to see that u has large tails, above the minimum decay 
estimate of (|9.1I1 . Then conservation of mass is true for t > r > 0. 

In order to get conservation of mass since the beginning we add a small perturbation with 
suitable tail to uq. After an easy argument this implies general conservation of mass and 
general minimum decay estimate. 

• General data with compact support. We attack the general case for data with 
compact support in the interval Ir = [—i?, R]. We propose to use convolution with a 
smooth kernel to obtain a smooth approximation uqs, that produces a solution to which 
the previous paragraph applies. By contraction we get conservation of mass, so that the 
solution must be global in time. 

Due to mass conservation and the smoothing effect the mass of the solution u cannot be 
contained in the original interval for large times. Then there is a time Ti such that half 
the mass is outside Ir. By the space monotonicity away from Ir, we can put a displaced 
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rearranged subsolution below, and then there are tails with at least the minimal decay rate 
for t > Ti. Now we use the monotonicity in time to derive the same conclusion also for all 
t >T > 0. 

• Finally, when uq is not compactly supported, rearranged or continous, we just approximate 
from below with compactly supported data and pass to the monotone limit. The rest of 
the argument follows. 

9.1 continuity and initial data 

We want to show that the limit solution is continuous as an orbit t i—)■ as stated in 

Theorem o At t = 0 this means that it takes the initial data in L^(M). We may use the 
fact that this holds for the approximate problems and then pass to the limit, using a Fatou 
argument plus conservation of mass. 

The continuity at t > 0 can be made into a stronger result, and in fact we obtain Lipschitz 
continuity of the orbit for all positive times. This is a consequence of the time monotonicity 
and the conservation of mass. Indeed, the first implies that for h > 0 

u{x,t + h) — u{x,t) < ((1 + — l)u{x,t) < Chu{x,t) 

which is uniform if t > r > 0 and h < ct. Since u is bounded, this implies a pointwise 
Lipschitz bound from above. The bound from below comes from conservation of mass. 

9.2 Upper limit solutions are very weak solutions 

This subsection extends to the singular case the result of papers [igiiniiso]. Since the ap¬ 
proximate solutions Us of Problem (12.ip are smooth, they satisfy the very weak formulation. 
The main difficulty in passing to the limit when e ^ 0 is the control of the possible growth 
of as |x| —>■ oo. But this depends on having a good lower bound for u, and such a 

bound is contained in the left-hand side of (j9.ip . It follows that the integrals involved in 
the passage to the limit are uniformly absolutely integrable. 

Now we can pass to the limit is the £ approximate equations written in very weak form to 
show that the limit solution is indeed a very weak solution. 

9.3 Upper limit solutions form a semigroup 

We consider the maps St : L^(]R) —>■ L+(M) that map any initial data uq to the solution of 
the equation at time t > 0, StUo = Since the approximate problems ()2.1I1 produce 

unique classical solutions, the semigroup property is true for them, = S'f(5f^uo). 

Passing to the limit we obtain the same property for our upper limit solutions arguing as 
follows: 

Fix no and ti > 0 and recall the monotone convergence Us(x, ti) —>■ n(x, ti) that takes place 
in Consider now vq = u{-,ti) as initial data for a new lap of the evolution. In order 

to get the upper limit solution St{vo) we take e' > 0 and solve the approximate problem with 
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vo,e' = vq + e' . If we now take s' smaller than e/2 we get an estimate of — ti))+||i 

in the following way. We first note that set of points K where u{x,ti) > e/2 has measure 
less than 2||n(x,ti)||i/e. Therefore 



Ue{x,ti))+ dx 


'K 


{vo{x) +e' 


Us{x,ti))+ dx < \K\e'. 


On the other hand, since Ue{x,ti) > e everywhere, and vq^^i{x) < e/2 + e' on R/ifT, on that 
set we have — Ue{x,ti) < 0. We conclude that 


s 

||(uo,£' - ?r£(x,ti))+||i < -||uo||i ■ 

By the ordered contraction property, this estimate remains true during the evolution, hence 

s' 

||(5’fuo,£' - Us{x,t + ti))+\\i < —||uo||i 

holds for all t > 0. Let now e' —)• 0 to get by the very definition of the upper limit solution 
that 

||(5tUo(x) -tte(x,t + fi))+||i < 0. 

This means that StVo{x) < Ue{x,t + ti)) for a.e. x G R. In other words, S't(S't^no) < St+tiUo- 
By the conservation of mass, both functions are the same, which proves the semigroup 
property St{Stj^uo) = St+tiUo- □ 


9.4 Upper limit solutions are mciximal solutions 

In the standard Laplacian case s = I with singular diffusion, it is known that in the range 
of exponents 0 < n < 1 where there is existence of nontrivial solutions, the limit solution 
is not the only possible solution defined in the whole line x G R. On the contrary, there 
are infinitely many other solutions with the same initial data determined by some “flux at 
infinity”, as described in |21l [28] . But the upper limit solution is the maximal element in 
that class, in fact it is the only one for which the total mass does not decrease in time. 

In order to repeat the proof of maximality in a short way, we only to consider a class of 
solutions that admits comparison with classical supersolutions (which will be the solutions 
Ue of the approximate problems). Let us call good solutions the elements of such a class. 

Theorem 9.2. Let Ug{x,t) denote a good solution of (jl.ll) for our choice of ^ defined in 
an interval 0 < f < T and having nonnegative initial data uq G L^(R), and let u{x,t) the 
upper limit solution with same initial data. Then, < U in R x (0, T). 

The further exploration of the existence of such solutions falls out of the scope of this 
paper for reasons of space. 
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10 Asymptotic behaviour of general solutions 

We want to prove Theorem 11.41 which means that we want to show that the Barenblatt 
solutions are asymptotic attractors of the solutions with general data restricted only by the 
running conditions on the initial data: uo{x) > 0, uq € L^(M). We address first the question 
of convergence in the norm, which is split into a number of cases. 

• We consider first rearranged initial data. The result is just a reformulation of the proof 
of construction of the Barenblatt solution in Section [6j The argument is well-known in 
the literature, see [38]. We just recall that the rescaled family UL{x,t) converges to Um in 
L^(M) at any time t > 0, and fix t = 1 to get 

lim \\uLix,l) - UMix,l)\\i ^ 0 
L—^oo 

We then undo the scaling, see formula (|5.5p . to find that 

lim \\u{x, ^ q 

L—>-oo 

Putting t = ^ 2 s-i-n ^ QQ obtain the result. 

• The main novelty lies in establishing the result for general data that are not rearranged. 
We have to use the trick introduced in paper [27] with Portilheiro for general data with 
compact support. We argue as follows: we fix ti 1, and assume the support of uq be 
included in [—i?, i?]. Define 

uiir) ■.= inf u{x,ti), uofr) := maxulx,ti). 

We easily verify that ui{r), U 2 {r) are nonnegative and radially symmetric functions, they 
are nonincreasing as functions of r for r > R, we have the immediate comparison 

ui{r) < u{x,ti) < U 2 {r ), 

where \x\ = r. We also have by the Aleksandrov reflection comparison 

U 2 {r) > ui{r) > U 2 {r + 2R) 

for all r > R. It is then easy to verify that the 1-d mass of U 2 {f’) — ui{r) is less than 
CRti which can be made very small. 

We restart the evolution at time ti and get radially symmetric solutions ui{x,t),U 2 {x,t) 
with initial data 'Ui(r),U 2 (r) resp., and we also have the original u{x,t + ti) (now displaced 
in time) that stays between them. The asymptotic behaviour says that ui{x,t) converges 
to the Barenblatt Umi, and ui{x,t) converges to the Barenblatt Um2- Moreover, the masses 
satisfy Mi < M < M 2 and M 2 — Mi < e. The asymptotic formula for convergence of u{-,t) 
to follows easily. 

• If no does not have compact support we use approximation and contraction. 

• The estimate in the U’ norms is just an interpolation between the convergence result for 

p = I just proved, and the bound of the form u{x, t) < C t~°‘ with a = l/(2s —1—n). □ 

This result allows to extend the uniqueness of the fundamental solutions as follows 
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Theorem 10.1. Every upper limit solution for positive times that that is self-similar and 
integrable is a Barenblatt solution UMix,t) for some M > 0. It can be obtained by rescaling 
from any integral and nonnegative initial data with mass M. 

Proof. Let ui be that solution and M its mass. Since it is self-similar we have 

- UM{{-,tl)\\l = ||'Ui(-,t2) - UM{{-,t2)\\l 

Fix now ti > 0 and let t 2 —oo. Applying Theorem 11.41 the right-hand side goes to zero. 
Hence, ui = Um- The second assertion is easier. □ 


11 Proof of the Shifting comparison lemma 

In this section we will give a complete proof of Theorem 14.21 Note that this result applies 
to the approximate problems that have nonsingular functions <1>. The passage to the limit 
allows to apply it to our upper limit solutions. 

11.1 Elliptic problem. Extended problem 

The implicit time discretization scheme mm directly connects the analysis of the parabolic 
equation (II.ip to solving a sequence of elliptic equations of the form 

(11.1) {-Ay/^v + B{v) = f{x) xGM, 

where a € (0,2) and / is an integrable function defined in M. We assume that the non¬ 
linearity is given by a function B : M_|_ ^ M_|_ which is smooth and monotone increasing 
with B{0) = 0 and B'{v) > 0. It is not essential to consider negative values for our main 
results, but the general theory can be done in that greater generality. We are using here 
nonnegative data and solutions. In the parabolic application B = see |38( 143]. Then 
we need to prove the following result 

Theorem 11.1. Let us consider two functions fi, f 2 G L^(M) with the following properties: 

(i) They are nonnegative, and rearranged around their points of maximum xi and X 2 resp., 
with xi < X 2 . BY this mean that fi{x — Xi) is rearranged. 

(a) We assume moreover that f 2 dx < /i dx 
(in) the mass is the same, f fi(x) dx = f f 2 (x) dx = M > 0. 

Besides, we assume that B is convex and defined on M+ with B{0) =0 and B'{0) > 0. 
Then, the following comparison inequalities hold 

/ X rx rx rx 

V 2 {x)dx< / vi{x)dx, / B{v 2 {x))dx< / B{vi{x))dx 

-OO j —OO J —oo J —oo 

for every x G M. 
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Proof, (a) It will be convenient to formulate the elliptic problem by using a proper extension 
problem, which is defined as the trace of a properly defined Dirichlet-Neumann problem as 
follows. If rc is a weak solution to the local problem 


(11.3) 


""Vw) =0 in Q+, 

1 d'Hj 

-lim -;—(x,y) + B{w(x,0)) = fix) for xGM 

y-^o+ oy 


where Q+ := Rx (0, +oo) is the upper half-plane and is the constant that is not important 
in what follows. See [TUjwhich is the main reference in the issue of this extension. We can 
define again a suitable meaning of weak solution in terms of this extended problem. The 
functional setting is perfectly explained in [13], Section 3.1. There is a solution to problem 
(lll.3p . Then the trace of w over R x {0}, TrR(r(;) = w{-,0) =: v is said a solution to problem 
(lll.ip . Using the change of variables 2 ; = Cy^ for a convenient constant c > 0, the problem 
can also be written as 


(11.4) 


^d‘^w d'^w 

(a:,0) = f{x) - B{w{x,f))) 


in Q+ 
for X G R. 


We use the extension formulation and write Wi{x,Q) = Vi{x) and Ui = B{vi), z = 1,2. 

(b) Note that the properties (i) of the ffs hold also for the solutions vi and V 2 . Conserva¬ 
tion of mass applies so that B{vi) and B{v 2 ) have the same mass M. Due to the properties 
of B the functions vi and V 2 are also integrable, though their masses are not controlled. By 
the properties of the extension to the upper half-plane, the functions W 2 {z, z) and wi{x, z) 
are also integrable in x for every fixed 2 ; > 0. 

(c) Using a similar strategy to the symmetrization proof in paper [43], we introduce the 
function 

(11.5) Z{s,z)= f {w 2 {t,z) - Wl{T,z))dT . 

J —00 


Then, it is clear that 


( 11 . 6 ) 


cz T Zxx — 0 


and 


(11.7) Z{—oo,z) = 0, Z{oo,z) = 0. 

A crucial point in our arguments below is played by the derivative of Z with respect to 2 ;. 
Due to the boundary conditions contained in (111.41) . we have 

(11.8) Z^{x,0)> f {B{w2{t, 0)) - B{wi{T,0))dT 

J —00 
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Recall that Wj(r, 0) = Observe also that the function 

y(x,0)= f B{w2{t,0)) - B{wit, 0)) dr 
J —oo 


has the same points of maximum or minimum and the same regions of monotonicity than 
Z{x, 0). 

(d) Then we argue as follows. Due to the maximum principle and the boundary conditions 
(lll.7p . a positive maximum of Z can be achieved only on the line {z = 0}. On the other 
hand, in the interval I = {x : xi < x < X2} we know that dxVi < 0 and dxV2 > 0 hence in 
this interval I 

Yxx = dxB{v2) - dxB{vi) > 0 

and the maximum of Y must lie outside of I, hence the same happens for Z. Suppose the 
maximum of Z happens at (xq, 0) with xq > X 2 - We must also have Zz{xo, 0) < 0 by Hopf’s 
maximum principle, and by (jll.81) . this leads to Y{xo, 0) < 0. But for x > xq 


y(x,o)-y(xo,o) = f 

s 

< f B'{v2{t,0)){v2{t) 

J Sn 


[B{v 2 {t)) - B{vi{T))]dT 
- vi{t)) dr. 


Here, we have used the convexity of B so that B' is an increasing real function and 


Biv2{T)) - B{vi{t)) < B'{v2{t)) (n2(r)) - ni(r)). 


After integration by parts in the expression for the increment of T, we get 
Y{x, 0) - y(xo, 0) < [B'{v 2 {t)){Z{t, 0) - Z(xo, 0))]" - 


B”{v 2 {t))V 2 ,x{t){Z{t,Q) - Z(xo,0))(iT. 


'XQ 


Since Z has a maximum at xq and B' is positive, the first term in the RHS is non-positive. 
As for the second, we have: B" > 0, V 2 ,x < 0, and Z{x,0) — Z{xo,0) < 0, hence the last 
term is also nonpositive. We conclnde that y(x,0) < y(xo,0) < 0 for all x > xq. This 
is a contradiction, becanse by the conservation of mass property at pins infinity we have 
y( 00 ,0) = 0. Therefore, there is no positive maximnm for Z on this side. 

(ii) Similar argnment on the other side, x < xi reversing the roles of vi and V 2 and the 
direction of integration, that starts now at + 00 . We conclude that Z{x,Q) < 0 everywhere. 

(hi) Once we have Z{x, 0) < 0 we also want to prove that Y (x, 0) < 0. We may use Lemma 
lll.2l below (a well-known resnlt), taking advantage of the convexity of B and choosing any 
convex, increasing function : [0, 00 ) —>■ [0, 00 ). This ends the proof of the comparison 
theorem in this case. □ 


Lemma 11.2. Let f,g € L^{Ll) be two rearranged functions on a ball D = R/j(0). Then 
f < g if and only if for every convex nondecreasing function $ : [0, 00 ) —)■ [0, cx)) with 
<h(0) = 0 we have 


(11.9) 


^{f{x))dx< / ^{g{x))dx. 
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This result still holds if R = oo and f,g £ with s' —>■ 0 as |x| —>■ oo. 

Remark. We have imposed severe conditions on the shape of uoi, ^* 02 ) and we have required 
<h to be concave. Are these conditions necessary? They are not for the same result with 
standard Laplacian instead of fractional Laplacian. 


12 A preview of logarithmic diffusion with s = 1/2 

This is a kind of exceptional case in the parameter diagram. The study is a bit different 
from the previous analysis, and much of our intuition comes from a similar special case that 
happens for standard diffusion s = 1, in dimension N = 2 with logarithmic diffusion n = 0. 
An explicit solution exists then and it was used in [35]. The formula is 

(12.1) U{x, t) = , X G 0 < f < T. 

(1 + \x\^y 

Note that U is positive only in the time interval 0 < t < T. As a very weak solution, it can 
be defined for all t > 0 and it vanishes identically for t > T. Another important property 
is the total mass decay with a constant rate 

1 p 

(12.2) — / U{x,t) dx =—Stt. 

dt Jr2 

But, as shown by a number of studies, this rate does not correspond the limit solutions, 
which are characterized by the rule dM{t)/dt = —dvr. A gap of non-uniqueness opens up 
and this is carefully described in [39], where related literature can be found 

12.1 Existence of a positive solution 

Luckily, in our one-dimensional case there also exists an explicit solution of the evolution 
equation 

(12.3) (-A)^/^(log?r) = 0 inMx(0,r) 

with a smooth initial condition u{x,0) G L^(M). It is given by the formula 

2A (T — t) 

(12.4) 17(x,t) = —^in R X (0,-Foo), 

with any A > 0, so it is indeed a whole family of solutions related by scaling. We see that 
U has the separate-variable type as in the previous example, and U{-,t) is in L^(R) for any 
0 < t < T. It is very peculiar that the solution becomes identically zero in finite time. This 
is the so-called hnite-time extinction phenomenon which is typical of some ranges of fast 
diffusion, see m for standard diffusion and [I9l|26] for fractional diffusion. 

In order to prove that (112. 4|) is a classical solution of the equation we write it in the form 
U{x,t) = (T — t)F{x) and we need to find an integrable profile F > 0 such that 

(12.5) (-A)^/V(x) = F(x). 
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It only remains to check that F{x) = 2/(1 + |xp) is a solution of this nonlinear elliptic 
equation. This is only a technical calculus result. 

Observe that the initial mass is in all cases (112. 4p 

r (It 

Mo = ||17o||i = 2r / ,^-^ = 27tT, 

J I + 

so that in the extinction time is given by T = ||?7o||i/27r. Accordingly, the mass decay rule 
is M'{t) = —27r, which is related to what happens for N = 2 and standard diffusion (see 
above). 

In the same way as in the previous analysis of the case s > 1/2, we can use this example 
and standard comparison to prove that for initial data uo{x) > c/(l + |xp) with c > 0, the 
limit solution is nontrivial, more precisely positive for 0 < t < T with T = c/2. We do not 
expect upper limit solutions to conserve mass since this does happen for the special case of 
standard diffusion s = 1, = 2, as described in [301 HU EH]- 


13 Comments, extensions, and open problems 

• We have given preference in the paper to the treatment of the more singular case of 
exponents n > 0. However, the logarithmic equation is also covered and the main results 
are proved to be true, as special case n = 0. This is shown not only at the formal level, but 
also technical details are given as needed. There is another variant of the limit case n = 0, 
the sign diffusion ut = Asign (u), that is related to the total variation flow and was treated 
in |6]. We are not covering the fractional version of this variant. 

• The actual behaviour of the solutions of logarithmic diffusion on the line with exponent 
s = 1/2 is a question to be investigated, and this will done in a separate work. For the 
non-singular equation dtu -|- (—A)®(1 -|- logu) = 0 the study was done in [20], see also [32] . 

• The existence of non-maximal solutions, that are not upper limit solutions and do not 
conserve mass, is an interesting open problem. For the standard Laplacian, it was solved 
in [211 EH] where it is proved that there are infinitely many solutions for every integrable 
initial data and they are determined by some flux conditions at infinity. A large related 
literature has developed, see e. g. [23l EH E8l E9] . 

• Equation (15.8p is a kind of fractional p-Laplacian equation for the mass function. It 
would be interesting to perform a study of its properties and applicability. 

• Elliptic problems and the Crandall-Liggett approach. A natural way to solve the evo¬ 
lution equation is by implicit discretization in time, a method that became basic in the 
early studies of the Porous Medium Equation, see the original paper CL71 or [351 Chapter 
10]. In the present situation, it means that we have to solve elliptic problems of the form 
(—A)^$(n) + u = f. Note that when <I>(u) = log(u), this equation takes the suggestive 
form (—A)®u + e" = f. Non-existence results are carefully described in [7] for the range of 
parameters 2s < n -|- 1. The existence theory in our parameter range is not difficult. We 
refrain from further details for reasons of space, but see [71 Section 9]. 
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• Dirichlet problem in bounded domains. In the case of the standard Laplacian, s = 1, no 
nontrivial solutions exist for the Dirichlet problem with zero boundary data. The question 
is open for our equations. 

• Problems with more general nonlinearities <I>. Symmetrization can be applied in some 
cases using the results of [H] to compare a general with the power cases we deal with. It 
can give the starting results on nontrivial existence. Non-existence results in that direction 
are explained in [7]. Of course, in the simple case where d>(u) = —cu“" (or <I>(u) = c log(M)) 
with some c > 0, the constant may be absorbed into the time variable and we need not 
make any changes to the theory. 

• Problems with other classes of initial data are worth studying. A simple example are 
the constant solutions do not belong to our class of integrable upper limit solutions. It is 
not difficult to construct solutions for initial data in L^(]R) with p > 1 using the smoothing 
effect and approximation by integrable data. We can consider in this way initial data that 
are not integrable with different decay rates at infinity. It is even possible to consider data 
that grow at infinity. 

• A different direction is establishing existence of upper limit solutions with nonnegative 
Radon measures as initial data. An example for that extension are the Barenblatt solutions 
that we have just constructed. Since the basic process for general measures is easy following 
the indications of [40l Section 4], we leave it to the interested reader. 
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